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Q '. ABSTRACT 

£: 

^5 Aims. We explore spectral energy distributions (SEDs), star formation, and dust extinction properties of galaxies in the Local Universe. 

, Methods. The AKARI All-Sky Survey provided the first bright point source catalog detected at 90 /jm. Starting from this catalog, we 

selected galaxies by matching AKARI sources with those in the IRAS PSCz. Next, we have measured total GALEX FUV and NUV 
flux densities by a photometry software we have specifically developed for this purpose. Then, we have matched this sample with 
SDSS and 2MASS galaxies. By this procedure, we obtained the basic sample which consists of 776 galaxies. After removing objects 

^ , with photometry contaminated by foreground sources (mainly in SDSS), we have defined the "secure sample" which contains 607 

galaxies. Using this galaxy sample, we have explored various properties of galaxies related to star formation and dust extinction. 
Results. The sample galaxies have redshifts S= 0.15, and their 90-^m luminosities range from 10 6 to 10 12 L a , with a peak at 10 10 Lq. 
The SEDs display a large variety, especially more than four orders of magnitude at M-FIR , but if we sort the sample with respect to 
90 /jm, their average SED shows a coherent trend: the more luminous at 90 p.m, the redder the global SED becomes. The M,--NUV - r 

■ color-magnitude relation of our sample does not show bimodality, and the distribution is centered on the green valley between the 

■ blue cloud and red sequence seen in optical surveys. We have established formulae to convert FIR luminosity from AKARI bands to 
' the total infrared (IR) luminosity Ltir. With these formulae, we calculated the star formation directly visible with FUV and hidden 

by dust. The luminosity related to star formation activity (Lsf) is dominated by Ltir even if we take into account the far-infrared 
• • ■ (FIR) emission from dust heated by old stars. At high star formation rate (SFR) (> 20 M G yr -1 ), the fraction of directly visible SFR, 

SFRpyv, decreases. We also estimated the FUV attenuation A FUV from FUV-to-total IR (TIR) luminosity ratio. We also examined the 
' Ltir/Lfuv-UV slope (FUV - NUV) relation. The majority of the sample has L TIR /L FUV ratios 5 to 10 times lower than expected from 

rS \ the local starburst relation, while some LIRGs and all the ULIRGs of this sample have higher Ltir/Lfuv ratios. We found that the 

attenuation indicator Ln R /L FUV is correlated to the stellar mass of galaxies, M t , but there is no correlation with specific SFR (SSFR), 
SFR/M,, and dust attenuation Ltir/Lfuv- 

Conclusions. Together, these results show that the AKARI FIS All-Sky Survey gives a representative sample of SF galaxies in the 
Local Universe. This sample will be a comprehensive standard of various properties of SF galaxies to be compared with, e.g., high-z 
SF galaxies. 
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1. Introduction 

Star formation history of galaxies is one of the most important 
and exciting topics in extragalactic astrophysics and observa- 
tional cosmology. Especially, exploring the "true" absolute value 
of the cosmic star formation rate (hereafter SFR) has been of a 
central importance. 

However, it remained difficult for a long time because of 
dust extinction. Even in the Local Universe, there aresome prob- 



lems to estimate the absolute value of SFR density because of 
different de pendence of various SFR e stimators on dust atten- 
uation (e.g.. lHopkins & Beacoml I2006L and references therein). 
Active star formation (SF) always comes along with dust pro- 
duction, because of various dust grain formati on processes re- 
lated to the final stage of stellar evolution (e.g..lDwek & Scalol 
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U98rtlDwe1dll99aiNozawa et all 120031: iTakeuchi et al.Ll2005d) 
Observationally, SFR of galaxies is measured by the ultravio- 
let (UV) luminosity from massive stars because of their short 
lifetime (~ 10 s yr) compared with the age of galaxies or the 
Universe. However, the UV photons are easily scattered and ab- 
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sorbed by dust grains. Hence the SFR of galaxies is always in- 
evitably affected by dust which is produced by their SF activ- 
ity. Since the absorbed energy is re-emitted at far-infrared (FIR) 
wavelengths, it is essential to observe galaxie s both at UV and 
FIR to have an unbiased view of their S F (e.gjBuat et al.Ll2005b 



ISeibert et al.ll2005HCortese et aUl2006tlTakeuchi et alll2005al) . 

To know the history of SFR in the Universe, we must know 
the starting reference value, i.e., the SFR density in the Local 
Universe: otherwise, we cannot define how much larger the SFR 
was in the past. However, since the volume of the Local Universe 
is limited by definition, an all-sky survey is the only viable way 
to refine our knowledge of local galaxies. On the "directly visible 
SF" side, the adven t of the UV surveyor-type satellite GALEX 
dMartin et al.ll2005l) has changed the situation of UV astronomy 
drastically. GALEX is performing an all-sky survey (All Sky 
Imaging Survey: AIS) at FUV (1530 A) and NUV (2300 A) wit h 
detection limits of 19.9 and 20.8 mag dMorrissev et all 120071) . 
as well as deep surveys in some selected regions. In a previous 
study, we have shown that GALEX AIS provides us with an un- 
precedent ed opportunity to explore the visible SF in the Local 
Universe dBuatet al.Ll2007al) . 

As for the "h idden" side of SF, the In frared Astronomical 
Satellite (IRAS; iNeugebauer et all 11984 has brought a vast 
amount of statistics of dusty galaxies in the Lo cal Universe 
by IR AS Point Source Catalog (PSC) (see, e.g., ISoifer et all 
1 19871) . Subsequently, FIR facilities with much higher sen s itivity 
have been launched, like ISO (e.g., iGenzel & Cesarskvl |2000t 
IVermaet al.ll2005l) and Spitzer (e.g. Jsoifer et alll2008l) 7and re- 
vealed the deep IR universe, but the latter two were observatories 
dedicated to pointed observations. 

In contrast to the latter two facilities, the Japanese IR satellite 
AKARI has performed variou s large-area surveys at IR wave- 
lengths dMurakami et alll2007l) after IRAS, especially including 
all-sky surveys at FIR and MIR. In particular, wit h the aid of the 
Far-Infrared Surveyor (FIS: iKawada et all [2007) onboard, var- 
ious IR surveys were performed by AKARI. AKARI FIS has 
four FIR wavebands centered on 65, 90, 140, and 160 /mi, and a 
FIR all-sky survey was completed by this instrument. Since the 
latter two bands are longer than 100 /mi which was the longest 
wavelength band of IRAS, the obtained sample of dusty galax- 
ies is less biased than the IRAS sample, i.e., thanks to the better 
sensitivity to cooler dust than IRAS, AKARI can detect galax- 
ies with dust emission with lower temperatures. Thus, AKARI 
is a very promising facility to bring new knowledge of galaxies 
with cold dust. In addition to these SF-related wavelengths, we 
need other various wavelength bands to examine physical prop- 
erties of galaxies, e.g., stellar mass [closely related to near-IR 
(NIR)], and intermediate stellar population (related to optical). 
For the former, we ha ve a set of all-sky data from 2-Micron All- 
Sky Survey (2MASS: ISkrutskie et ail 120061) . and for the latter, 
the SDSS final data (DR7) are publicly available^, even if SDSS 
is not an all-sky survey. 

In this work, we constructed a multiband galaxy catalog 
based on AKARI All-Sky Survey 90-/mi selected source s asso- 
ciated with IRAS PSCz galaxies dSaunders et all |2000|) . Then, 
we measured GALEX FUV and NUV flux densities, and asso- 
ciated SDSS and 2MASS photometries. For this initial study, 
we have only selected "secure" galaxies with good photometric 
measurements for most of the bands. We present the sample con- 
struction in Section [2] We describe basic properties of galaxies 
in the sample in Section|3] Results on the SF and dust attenuation 
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Fig. 1. Comparison between AKARI FIS and IRAS PSCz flux 
densities. Upper-left, upper-right, and lower-left panels present 
comparisons of IRAS 60 //m with AKARI 90 /mi, 65 //m, and 
140 /jm flux densities of the AKARI-IRAS correlated sample. 
The vertical dotted lines in these panels represent the flux density 
limit of IRAS PSCz. Lower-right panel shows a comparison of 
IRAS 100 /mi with AKARI 90 /mi flux densities. The horizontal 
dotted line in Upper-left panel represents the formally expected 
detection limit of FIS 90 /mi. 



are presented in Section|4] Section[5]is devoted to our summary 
and conclusions. 

Throughout the paper we will assume Qmo = 0.3, Q.ao = 0.7 
and Ho = 70 km s _1 Mpc -1 . The luminosities are defined as vL v 
and expressed in solar units assuming L = 3.83 x 10 33 erg s _1 . 



2. Sample construction 

2.1. Construction of the parent FIR sample from AKARI and 
IRAS PSCz 

2.1 .1 . AKARI FIS All-Sky Survey 

The primary purpose of the AKARI mission is to provide 
second-generation infrared (IR) catalogs with better spatial res- 
olution and wider spectral coverage than the IRAS catalog. 
AKARI is equipped with a cryogenically cooled telescope 
of 68.5 cm aperture diameter a nd two scienti f ic inst ruments, 
the Far-Infrared Surve yor (FIS; Kawadaetal., 2007) and the 
Infrared Camera (IRC; lOnaka et al.Ll2007T) . Among various as- 
tronomical observations performed by AKARI, as we have men- 
tioned in Introduction, an all sky survey with FIS has been car- 
ried out (AKARI All-Sky Survey). Since FIS is an instrument 
dedicated to FIR A = 50-180 /mi, all the AKARI FIS bands 
are in the FIR wavelengths: N60 (65 /m i), WIPES (90 /mi) , 
WIDE-L (140 /an), and N160 (160 /mi) dKawada et all 120071) . 
Hereafter, we note as $65, Sgo, S ho and S i6o the flux densities 
in these bands. Especially, since FIS has sensitivity at longer 
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wavelengths than IRAS, a new classification scheme is needed 
if we try to select a certain class of objects. Such a scheme is 
not merely an empirical technique but also will provide us with 
a new understanding of objects with cool dust which were diffi- 
cult to detect by IRAS bands. 

We use the AKARI FIS Bright Source Catalogue (BSC), 
the first primary catalog from the AKARI All-Sky Survey. Data 
from the version of this catalog are used in this work. 
AKARI BSC is supposed to have a uniform detection limit, cor- 
responding to per scan sensitivity, over the entire sky, except for 
very bright sky parts where different data acquisition mode had 
to be applied. A summ ary of the All-Sky Survey is presented in 
lYamamura eTail d2009l) . The AKARI FIS BSC provides data for 
6431 1 sources. For each detected source, AKARI source identi- 
fier, equatorial coordinates of the source position and flux densi- 
ties in the four FIR bands are given. Errors are not estimated for 
each individual source, but instead they are in total estimated to 
be 35 %, 30 %, 60 %, and 60 % at N60, WI PES, WIDE-L, and 
N160, respectively ( lYamamura et all 120081) . AKARI IRC per- 
formed another all sky survey, but the data are not fully reduced 
at the time we have been preparing this paper. Hence we focus 
on the FIS data only. 

2.1.2. Matching with IRAS PSCz 

The AKARI BSC contains many Galacti c sources, l i ke AG Bs, 
Hn regions, planetary nebulae, etc. (e.g., iPollo et all |201Q|) . In 
order to construct a reliable catalog of galaxies, we should pick 
up FIS sources confirmed as galaxies. For this purpose, we have 
perfor med a cross identificat ion of BSC sources with the IRAS 
PSCz dSaunders et all |2000t) . The PSCz is a redshift survey of 
galaxies selected at IRAS 60 fim with a flux density limit of 
Seo > 0.6 [Jyfl The PSCz contains ~ 16000 galaxies. We have 
put a limit on recession velocity v > 1000 kms" 1 so that we can 
avoid the effect of the peculiar velocity of galaxies. Then, we 
have matched the AKARI BSC sources with PSCz galaxies with 
a search radius of 36 arcsec, which corresponds to the position 
uncertainty of IRAS PSC. The number of matched sources was 
5890. To examine the effect of the choice of search radius, we 
changed the criterion from 20-60 arcsec. This change of radius 
does not have a significant impact on the resulting catalog (^ 
5 %). Hence, we use the catalog with a search radius of 36 arcsec 
in the following analysis. 

For further analysis, we make a cross identification with 
SDSS galaxies (see Sec. |2.31 >. Then, we have to restrict our data 
only to the area covered by SDSS DR7 which is 8378.015 deg 2 . 
This restriction to this region of the sky has one advantage: the 
Galactic diffuse FIR emission is strong in some areas of the sky. 
In such regions, measured FIR flux densities of point sources are 
contaminated by the Galactic emission and not very accurate. 
Since the SDSS region is selected so that the Galactic extinction 
is small, the selected area automatically avoids such FIR-bright 
regions. Then, by this selection, our sample automatically ex- 
cludes strongly contaminated sources. The resulting parent cat- 
alog contains 1186 galaxies. 

2.1.3. Flux density comparison between AKARI and IRAS 

We compared AKARI and IRAS flux densities to examine our 
sample selection. The correlation is shown in FigureQ] The hori- 



2 Because of this step, we should note that we do not make a maximal 
use of the advantage of the long wavelength bands of AKARI FIS, since 
the sample is limited by IRAS bands (A < 100 pm). 



zontal dotted lines in upper-left, upper-right, and lower-left pan- 
els represent the flux density limit of IRAS PSCz. AKARI BSC 
sources are selected at WIDE-S, i.e., 90 fim. It is important to 
see which selection controls the sample selection. As seen in the 
upper-left panel, the IRAS PSCz limit and AKARI limit are both 
well-defined, and neither of them strongly restricts the sample. 
The effective 90 fim flux density limit of our parent sample is 
~ 0.8 Jy. 

2.2. GALEX photometry 

GALEX AIS now observed 25000 deg 2 at FUV and NUV. The 
latest version of the public imaging is GR4/GR53 We have 
measured the FUV and NUV photometry of the parent AKARI 
galaxies as follows: 

1 . Cut out a 30' x 30' square subimage from GALEX AIS im- 
ages around each AKARI galaxy. 

2. Select a subimage with the largest exposure time when mul- 
tiple observations were available. 

3. Measure FUV and NUV flux densities. The NUV observa- 
tion is taken as our reference. 

Since the sky coverage of GALEX AIS is not complete, in some 
cases we do not have a proper GALEX image for an AKARI 
galaxy. In such a case we omit the galaxy because we do not 
have any UV information. 

Almost all of the sources are resolved by GALEX. They 
are thus very often separated into small bright patchy regions, 
and the GALEX pipeline misidentifies these fragments as indi- 
vidual objects. This is referred to as shredding. We must deal 
with the shredding to obtain sensible flux density measurements 
for nearby extended galaxies. For this purpose, we have used 
an IDL software package developed by ourselves. This software 
performs aperture photometry in the NUV sub-image using a 
set of elliptical apertures. Total flux density is calculated within 
the aperture corresponding to the convergence of the growth 
curve. The sky background is measured by combining several 
individual regions around the source. NUV and FUV flux den- 
sities are_TOrrec_ted forGalactic extinction using the Schlegel 
map dSchlegel etal| [T998b and the Galactic extinction curve of 
ICardelli et alHim% . A detailed description of the p hotometry 
process can be found in llglesias-Paramo et all d2006l) . 

This software was already used for our previous IRAS- 
GALEX based studies dBuat et all l2007al) . and its performance 
is carefully checked and established. During the procedure, we 
also excluded all the sources contaminated by stars or too close 
to be disentangled by the photometry software. By this proce- 
dure, the number of remaining galaxies is 776. 



2.3. Matching with SDSS and 2MASS galaxies 

We further matched the AKARI-IRAS PSCz-GALEX sample 
(776 galaxies, hereafter abbreviated as the AKARI-GALEX 
sample) with SDSS and 2MASS. 

The AKARI-GALEX sample was matched with the 2MASS 
All-Sky Extended Source (XSC) catalog in order to obtain the 
NIR (/, H, and Ks) flux densities. A matching radius of 20" was 
initially adopted for the cross correlation. All but 7 out of the 776 
galaxies in the AKARI-GALEX sample showed 2MASS coun- 
terparts at distances closer than the matching radius. In case of 
multiple candidates, we always selected the brightest one at Ks 
band as the most plausible one. A further cross correlation with 
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Fig. 2. Number counts of our AKARI multiband sample at Fig. 3. Redshift distribution of the sample. The distribution is 

AKARI FIS four bands. Triangles, diamonds, squares, and normalized so that we obtain the number of galaxies per unit 

crosses represent AKARI N60, WIDE-S, WIDE-L, and N160 solid angle if we integrate it over redshifts. 
galaxy counts, respectively. 
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a matching radius of 30" was attempted for these seven galax- 
ies with no 2MASS counterpart within the 20" radius. However, 
again, no 2MASS counterparts were found, and we adopted the 
standard 2MASS XSC upper limits for these galaxies at J, H, 
and Ks bands. The standard 2MASS limiting magni tudes at J, H , 
and K s are 14.7, 13.9 and 13.1 mag, respectively dJarrett et al.L 
120001) . 

We also matched the AKARI-GALEX sample to SDSS us- 
ing the GALEX coordinates for the AKARI-GALEX objects, 
and a search radius of 15". In the case of SDSS, large resolved 
galaxies such as those we are dealing with here are shredded in 
multiple detections during the deblending step of the pipeline. 
We used the closest SDSS match to the AKARI object to obtain 
the SDSS photometry of the parent object, namely the object 
detected by the SDSS pipeline before deblending. We inspected 
all sources to check that the actual flux measured for the parent 
object is not contaminated by nearby bright stars, artifacts etc. 

For most of the sample galaxies, SDSS galaxies were as- 
sociated. However, stars superposed on the SDSS galaxy im- 
ages often hamper accurate photometry. It requires a very care- 
ful masking of the SDSS image, and in the current analysis we 
simply omitted such galaxies. After this selection, 607 galax- 
ies remained. We call this sample "the final sample". Almost all 
galaxies have very secure photometric data at UV, optical, NIR, 
MIR, and FIR, as well as the redshift information. 



3. Basic Properties of the Sample Galaxies 

3.1. Number counts 

Figure [2] presents the cumulative AKARI flux density distribu- 
tion (number counts) of our final sample. We see that the flux 
density limit of the sample at 90 /im is ~ 1 Jy, shallower than 
the original limit of ~ 0.8 Jy, because of additional constraints 
to have optical-NIR counterparts. 
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Fig. 4. Luminosity distribution of the sample at 90 fim. 



3.2. Redshift and luminosity distributions 

The redshift distribution of our final sample is shown in Figure[3] 
Since AKARI FIS BSC is rather shallow, most of the sam- 
ple locate at low redshifts z 5= 0.05. This is quite consistent 
with predictions of va ri ous IR galaxy evoluti on models (e.g., 
ITakeuchi et ail I200lallb1: ICharv & Elbazll200Tl) . 

Figure [4] presents the 90-/mi luminosity distribution of the 
sample. Since this is a raw luminosity distribution of IR galaxies, 
it decreases toward the faint end. The 90-/im luminosity ranges 
from 10 6 Lq-10 12 L q . Only a few galaxies are classified as ul- 
traluminous IR galaxies (ULIRGs). The peak of the luminosity 
distribution is around 10 10 L Q , which is more luminous than th e 
knee of the 60-/Ytn luminosity function ( ITakeuchi et all l2~003al) . 
Since our selection procedure is multifold and complicated, esti- 
mating a reliable luminosity function is not straightforward. We 
will try this task in future works. 
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Fig. 5. Luminosity distribution of the sample at several wave- 
lengths. Triangles, diamonds, squares, and crosses represent 
galaxies AKARI WIDE-S, GALEX FUV, SDSS r, and 2MASS 
/if -band, respectively. 



We can compare the distribution of galaxies on the redshift- 
luminosity plane to see the effect of the shape of luminosity 
function. We show the z-vL v relation in Figure [5] Some rep- 
resentative wavelengths are shown: AKARI WIDE-S, GALEX 
FUV, SDSS r, and 2MASS Ks. We see that the luminosity at 
90 /im increases monotonically toward higher redshifts, while 
the optical and near-IR (NIR) luminosities of the same galax- 
ies saturate at certain values. This is related to t he difference in 
the shape of their luminosity function s (see, e.g., Takeuchi et al.L 
l2005ah Uglesias-Paramo et all , 120061) : i.e., the luminous end of 
the function exponentially declines at optical-NIR and UV, while 
it shows a power-law decline at FIR. This is also related to the 
fact that the more luminous galaxies are more strongly extin- 
guished by dust. We will revisit this issue in a much more direct 
way in Section l4~4l 



3.3. SEDs of the sample 

Since we have the monochromatic luminosities from FUV to 
FIR, we construct the spectral energy distributions (SEDs) of our 
sample. We show all the SEDs of the sample in Figure[6] Large 
empty squares represent the AKARI FIS measurements, while 
the small squares represent all the other data in Figure[6] We ob- 
serve a very large variety of SEDs among the sample galaxies. 

To see a global trend of the SEDs more clearly, we sorted 
the sample by their 90-/im luminosities and subdivided the sam- 
ple into six logarithmic bins from 10 6 L Q to 10 12 L (with a 
bin width AlogLgo = 1)- By taking an average and median, we 
constructed "average SEDs" as a function of 90-fim luminosity, 
shown in Figure Q Because of a large dispersion and asymmet- 
ric distribution of the SEDs in each bin, sometimes the average 
and median SEDs do not agree very well. Even so, now the trend 
is more clearly seen: low 90-pm luminosity galaxies have cooler 
dust emission and bluer UV-optical continuum, while high 90- 
fim luminosity galaxies have hotter dust emission and redder 
UV-optical continuum. 

The 90-fim luminosity dependence of the dust emission tem- 
perature is more clearly seen if we plot a flux density ratio 
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Fig. 6. Spectral energy distributions (SEDs) of the whole sample. 
Large empty squares represent the AKARI measurements, while 
the small squares are data taken from GALEX FUV, GALEX 
NUV, SDSS u, g, r, i, z, 2MASS J, H, Ks, and IRAS 12, 25, 60, 
and 100 fim from left to right. Dotted lines connect data points 
of each object to guide the eye. 



S 65 /S 90 as a function of L90. This is shown in Figure [8] Again, 
the agreement between average and median values is not excel- 
lent, we see a clear monotonically increasing trend of the flux 
density ratio along with L90. 

The lowest luminosity galaxies have the largest uncertainty 
mainly because of poor statistics. This problem will be solved by 
constructing deeper and larger sample, possibly by the next gen- 
eration AKARI FIS catalog. The most 90-/im luminous galaxies 
seem to have an upturn at UV. This may be because of an AGN 
component in these galaxies. We will examine this issue in our 
future work. In both Figures|6]and|7] it is clear that AKARI FIS 
measurements at N60 and WIDE-S do not agree with IRAS ones. 
This discrepancy is mostly due to the better angular resolution 
of AKARI FIS compared with IRAS. Thus, the measured values 
are much less contaminated by Galactic cirrus emission. Also 
because of better angular resolution, the source confusion effect 
is much smaller than for IRAS. As both cirrus and source confu- 
sion effects cause flux boosting for IR galaxy n umber counts, the 
IRAS flux densities coul d be overestimated dTakeuchi & Ishiil 
120041) . iJeong et al.l(l2007h also examined the same problem with 
an early AKARI sample and concluded that the difference be- 
tween IRAS and AKARI flux densities are due to the confusion 
effects. We see a strange bump at 12 fim for luminous galaxies. 
This may not be a physical feature in the SEDs of our sample 
but because of a poor measurement at this band by IRAS. In 
the redshift-luminosity diagram at 12 /im (see Fig. IA. Il l, we see 
that many galaxies locate on the limiting luminosity line. This 
means that, even if they are classified as measurements, actually 
they are upper limits of the flux density. We plan to study this 
point further by AKARI IRC all-sky survey in the future. 
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Fig. 7. Averaged spectral energy distributions of our sample as 
a function of 90 /jm luminosity. Luminosity bins used here are 
the same as the bins in Figure [4] The central luminosities of the 
bins are 10 7 , 10 8 , 10 9 , 10 10 , 10 11 , and 10 12 L , respectively, and 
the obtained average SEDs are represented with different size 
of the symbols: from the smallest corresponding to 10 7 L Q to 
the largest corresponding to 10 12 L Q . Filled and empty squares 
are averaged SEDs, while filled and empty triangles are median 
SEDs. 
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Fig. 8. Dependence of the flux density ratio S65/S90 on aver- 
age 90-fim luminosity. Same as Fig. [7] the squares are estimated 
from averaged SEDs, while triangles are from median SEDs. 



3.4. The NUV- r distribution 

The NUV - r restframe color is very ef ficient to separate 
the galaxies into b lue and red populations dSalim et al.L 120071: 
iMartin et all 120071) . At high-z, IR selected galaxies observed 
with Spitzer were found to populate the so-called green valley 




-23 -22 -21 -20 -19 
? — band absolute magnitude M r 



Fig. 9. The absolute r-magnitude-NUV - r color distribution. 
Upper and lower thick solid lines r epresent the averag e red se- 
quence and blue cloud taken from ISalim et all d2007l) with its 
envelope indicated by dotted lines. 



between the blu e cloud and red sequence, first mentioned by 
iBell et alJ d2005l) and further identified in t he COSMOS fields 
dKartaltepe etaUl2009llv"ergani et al.Ll2009t) . 

However, in contrast, no systematic study of IR selected 
galaxies has been performed in the nearby universe until now. 
Figure [9] shows the distribution of the NUV - r color for our 
sample galaxies against the absolute magnitude at r-band, M r . 
The distribution in NUV - r is found to be large and mono- 
lithic in contrast to the bimodal distr i bution found in optical sur- 
veys. From Figure 1 of ISalim et ail d2007l) (based on an SDSS- 
GALEX selected sample), we approximate the red sequence on 
this diagram with a linear relation 



NUV - r = -0.143M, + 2.57 , 

with ±0.5 envelopes, and the the blue cloud with 

NUV - r = -0.429M, - 6.29 , 



(1) 



(2) 



with +0.7. Clearly, we see that our galaxies lie in the blue 
cloud and populate the green valley to produce a continuous 
distribution, and only a few objects are located in the red se- 
quence. Therefore, the color distribution of FIR-selected galax- 
ies is different than the one obtained from an optical selection, 
as they preferentially populate the green valley between the red 
sequence and blue cloud. An interpretation of this behavior in 
terms of dust attenuation, SF history or AGN activity will be 
developed in a future paper. 



4. Star Formation and Dust Extinction of the 
90 //m-Selected Galaxies 

4. 1 . Total dust luminosity from AKARI FIS bands 

To calculate SF- and attenuation-related physical properties 
of the galaxies, a total IR (TIR) luminosity is required. 
Various estimators of Ltir were proposed by previous au- 

JPale & Heioull2002l: 



ire p 

thors dHelou et all|l988tlf>ale et aUl2001. 
ISanders & Mirabeli I1996I) . iTakeuchi et al 



(2005b) has shown 
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Fig. 10. Relation between the AKARI FIR luminosity and total IR (TIR) luminosity of the sample. Left panel is for the AKARI 
FIR luminosity defined by eq. ©, while right panel is for the AKARI FIR luminosity from 2 bands, defined by eq. ([8]). Black solid 
lines are the best least-square fits to the data, and the dotted curves represent the 95 % confidence levels of the lines. Triangles are 
galaxies with IRAS 100 /mi quality flag larger than 3, i.e., those with insecure flux density measurement. 



that we c an safely estimate the TIR luminosity by a formula pro- 
posed bv lSanders & Mirabell (1 19961) : 



Ltir = 4.93xl0" 22 [13.48L v (12yum) + 5.16L v (24^m) 
+2.58L v (60/mi) + L v (100/mi)] [L Q ]. 



(3) 



Since our sample contains all IRAS band flux densities, we can 
use the estimator of Sanders & Mirabel to examine the AKARI 
FIR flux. In this subsection, we make an attempt to establish a re- 
liable formula to convert a FIR luminosity measured by AKARI 
to the TIR luminosity. 

Since AKARI FIS has continuous bands from ~ 50 fim to 
~ 160 /im, we can easily define FIR flux simply by adding 
the flux densit i es mu ltiplied with their bandwidths (in [Hz]). 
iHirashita etaTI d2008l) has defined the AKARI FIR luminosity 
as 

Lakari = Av{N60)L v (65 fim) + Av (WIDE-S) L y (90 /urn) 

+Av (WIDE-L) Ly( 140 fim) , (4) 

where 

Av(N60) = 1.58 x 10 12 [Hz] (5) 
Av (WIDE-S) = 1.47 x 10 12 [Hz] (6) 
Av (WIDE-L) = 0.831 x 10 12 [Hz] . (7) 

However, since the sensitivity of AKARI N60 is not as good 
as other two wide bands, if we can use FIR luminosity defined 
only by WIDE-S and WIDE-L , it will be useful because we can 
have larger number of galaxies. We define LavTrt by omitting 
the term of N60 as 

t 2band 



,, = Av (WIDE-S) L v (90 fim) + 
Av (WIDE-L) L v ( 140 fim) . 



(8) 



Figure [10] presents correlations between Lakari and 
Ltir. Left panel sho ws the correlation between Lakari of 
IHirashita et ail (120081) and Ljir, while right panel shows the one 
between L A 1 £ A ^ U and Ltir. The fitting results are as follows: 

log Ltir = 0.940 log Lakari + 0.914 , (9) 
r = 0.987 , (10) 



and 

log Ltir = 0.964 log L 2b ^ A d RI + 0.814, 
r = 0.989 , 



(11) 
(12) 



where r is the correlation coefficient. The solid lines in Figure [TU1 
depict these equations. The envelopes delineated by dotted lines 
in Figure [TOlrepresent the 95 % confidence intervals. 

In some cases the quality of IRAS 100 fim flux density mea- 
surement is poor; we examine the impact of these objects in 
Figure[l0] Triangles represent galaxies with IRAS 100 /im qual- 
ity flag larger than 3, i.e., the ones for which the measurement 
was difficult and not secure. However, even so, clearly they lie 
well within the distribution of the whole sample and do not have 
a significant impact on the statistical analysis. 

Since these plots present luminosity-luminosity correla- 
tions, the very tight correlations are not extremely surprising. 
However, uncertainties of these equations are within a factor of 
~ 2-3, which means that the estimation works very well even 
if we do not have MIR measurements. We also stress that the 
correlation is even better if we only use AKARI wide bands. 

Hereafter, we use the total IR luminosity Ltir obtained from 
AKARI N60, WIDE-S, and WIDE-L [equation ©] for all the 
following analysis. Namely, when we mention Ltir, it is always 
AKARI-based total IR luminosity. 

4.2. Luminosity from current star formation 

Now we can compare the FUV (Lfuv) and total IR (TIR) (Ltir) 
luminosities for our final sample. This is shown in Figure QT| 
It is striking that the luminosity is dominated by Ltir f° r the 
vast majority of our sample, even though we take into account 
that they are FIR-selected. Also, it is worth mentioning that the 
luminosity at FUV does not exceed 10 11 L . In contrast, L T ir 
can be much higher. 

By combining Lfuv and Ltir, we can define the luminosity 
contribution from newly formed stars, Lsf: 



L SF = Lfuv + (1 - rj)L 



TIR 



(13) 
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Fig. 11. Comparison between TIR and FUV luminosities. The 
diagonal dotted line represents the case if Ltir equals Lfuv- 



where rj is the fraction of the dust emission due to the heating 
of grains by old stars which is not directly relat ed to the recent 
SFR. We adopt a value of 30 % for this fraction dHirashita et all 
12003b . The contribution of Lfuv and Ltir is shown in Figure [12l 
Naturally, the contribution of (1 - ^Ltir dominates Lsf- In con- 
trast, the contribution of Lfuv has a very large scatter, and the 
correlation is very poor. Hence, it is not surprising that it is very 
difficult to estimate the total energy emitted by newly formed 
stars only from FUV information, even in an average sense. 
The UV contribution becomes significant at lower luminosities 
Lsf < 10 10 L , as seen in the left panel of Figure [T2l 



4.3. Star formation rate of galaxies in our sample 

Here, we interpret the data in terms of SFR. Assuming a con- 
stant SFR over 10 8 yr, and Salpeter initial mass function (IMF) 
tiSalpeteri[l95l mass range: 0.1-100 M )), we have the relation 
between the SFR and L F uv 



logSFRpuv = log Lfuv - 9.51 



For the IR, to transform the dust emission to the SFR, we assume 
that all the stellar light is absorbed by dust. Then, we obtain the 
following formula under the same assumption for both the SFR 
history and the IMF as those of the FUV, 



4.4. Dust attenuation 

AAA . Dust attenuation of the sample 

Galaxies selected in IR are expected to have a quite large dust 
attenuation. Here we examine the extinction properties of galax- 
ies in the sample. For this study, a good observational indicator 
of dust attenuation is required. The Ltir/Lfuv ratio is widely 
recognized to be a robust measure of dust attenuation. This ra- 
tio was found to increase with the star f ormation lum i nosity 
in a similar way from z — to z — 0.7 (Martin e t all 120051 
iBuat et all l2007bh IZheng et ail 120061). Dust a ttenuation can be 
derived using the formula of IBuat etail d2005l) : 



Apuv = -0.0333 log 



+0.3522 log 



+ 1.1960 log 



i+UV 

Ltir \ 
Lfuv/ 
Ltir 



+ 0.4967 [mag] 



(17) 



Figure[l4]presents the variation of Ltir /Lfuv as a function of the 
star formation luminosity. A clear increase of Ltir /Lfuv with 
Lsf is seen. The solid li ne is the mean tren d of local galaxies 
found in our IRAS study dBuat et all 12007 al) : 



log 



= 0.64 log L SF - 5.5 . 



(18) 



This line was estimated by a weighting of L T ir/L F uv with 
1/Vmax to eliminate the flux selection effect, i.e., the effect that 
more luminous galaxies can be more easily detected. In contrast, 
we plot raw values of Ltir/Lfuv- However, though the steep rise 
of the distribution of Ltir/Lfuv with Lsf is partially because of 
this selection effect, the trend i s still well consistent with our 
IRAS study dBuat et all l2007al) . and the conclusion would re- 
main valid. 



4.4.2. Dust attenuation versus the slope of the UV continuum 

The slope of the UV continuum is commonly used as a proxy 
to estimate dust attenuation and Ltir/Lfu v when IR data are 
not availabl e or considered as unrel i able ( Daddi et ail 120071: 



(14) IReddv et all 120081: iReddv & Steidel. |200J). This method is 



based on a calibration pe rformed on local starburst galaxies 
dMeurer et all 119951 fl999l) . GALEX observations brought large 
amount of data in this field since the slope of the UV contin- 
uum can be safely deduced from the FUV - NUV color. Studies 
based on GALEX data have shown that the local starburst cal- 



log SFR dust = log Ltir - 9.75 + log(l - rj) 



(15) 



Here, again, t] is the fraction of the dust emission by old stars. 
Thus, the total SFR is simply 



SFR = SFRpuv + SFR dust 



(16) 



The obtained SFR is shown as a function of the fraction of 
the contribution of SFRpuv in Figure Qj] Reflecting the large 
scatter of Lfuv /Lsf, the scatter of SFRpuv /SFR is very large at 
SFR < 20 M yr~'. However, quite sharply, no galaxies have a 
large contribution of SFRuv at SFR > 20 M yr~'. The vertical 
dotted line shows this "threshold" SFR in Figure [T3l 



(e.g., Dale et al. 


2005tlCortese et all 


2006tlBoissier et alll2007t 


iGil de Paz et al. 


20071 Uohnson et al. 


, 2007): i.e., these galaxies 



exhibit lower dust attenuation (Ltir/Lfuv ) than expected from 
their FUV -NUV color. IR selected galaxies also depart from the 
starburst law but in the opposite way: they ar e found to exhibit 
a Ltir/Lfuv ratio mu ch larger than expected dBuat et all 120051 : 
iGoldader et all 120021) . Here, we re-inves tigate this issue m aking 
use of a much larger sample than used in lBuat et all d2005l) . 

Figure[T5lshows the FUV-NUV color against Ltir/Lfuv (of- 
ten called the IR-excess (IRX)-/? relation) for our sample galax- 
ies together wit h the local starburst relation. The local relation 
was taken from iMeurer et al. 1 (11999b and converted into the re- 
lation between the FUV-NUV color and Ltir/Lfuv by eqs. (1) 
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Fig. 12. Contribution of Ljir and Lpuv to the total star formation (SF) luminosity Lsf, which is the luminosity produced by newly 
forming stars defined by eq. (fTJl i. In the left panel, rj is the fraction of IR emission produced by dust heated by old stars which is not 
related to the current SF (see main text). 
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Fig. 13. Contribution of the FUV-estimated (or "directly vis- 
ible") SFR, SFR UV , to the total SFR as a function of total 
SFR. Vertical dotted line represents the effective boundary above 
which almost all energy produced by newly forming stars is 
emitted at IR. 

and (2) of lKong et al.l (12001 : 

^TIR _ jq[1.92(FUV-NUV)+0.4] _ q 95 nm 
^FUV 

Most of the galaxies lie below the starburst relation by a fac- 
tor of 5 to 10, but some locate above the relation. A significant 
fraction of galaxies above the local starburst line are luminous 
IR galaxies (LIRGs and ULIRGs). Especially, all the ULIRGs 
have larger Ltir/Lfuv ratios than expect ed from the relati on. 
This general trend is quite consistent with iBuat et alj d2005l) . It 
may be worth mentioning that LIRGs are consistent with the lo- 
cal starburst line. These trends are more clearly represented by 
functional fits in Right panel of Figure [15] We used a functional 



Fig. 14. Distribution of Ljir-^fuv ratio as a function of the 
star formation luminosity L c F . Solid line is th e relation found 
by IRAS-GALEX analysis feuat et all l2007al) . represented by 
eq. ([H. 

form of eq. ( fT9] l for "normal" galaxies with L^q < 10 11 L , 

^ T1R _ 1 Q[fl(FUV-NUV)+i] _ c ^0) 

^fuv 

with letting three parameters, a, b and c free. The fitted relation 
is well below the lo cal starburst relation originally claimed by 
iMeurer et alj ([1999). For more luminous galaxies, eq. ( f20b does 
not give a reasonable fit, hence we adopted a simple linear fit. 
Even if LIRGs are roughly consistent with the local starburst re- 
lation, the fitted line is much flatter than that, because of the exis- 
tence of galaxies well above the local starburst line with blue UV 
color. It is difficult to conclude in the case of ULIRGs because 
of poor statistics, but the fitted linear relation seems to locate 
above the local starburst relation. We have seen that galax- 
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ies with relatively quiescent galaxies (L(, p < 10 11 L ) in ou r 
sample locate well below the relation of iMeurer et al.l (1 19991) . 
Similar trend was reported for optically selected galaxie s (e.g. , 
iBoissier et all. 120071: ICortese et all 120061) IBoissier et all d2007l) 
investigated this relation for normal galaxies selected at opti- 
cal wavelength with the same functional form with this study 
[eq. (f20l >1 They found a much flatter relation than Meurer's orig- 
inal relation, which is pr esented in Right pane l of FigurePT"5lwith 
dot-dot-dot-dashed line. ICortese et al.l (120061) reported a similar 
result. What is newly clarified in this study is that the same trend 
is seen for purely FIR-selected sample of SF galaxies. In addi- 
tion, Boissier et al.'s relation does not represent the relation for 
our quiescent galaxy sample. This is naturally understood that 
our sample includes strongly reddened galaxies by dust which 
can never be picked up by UV-selection. 



4.4.3. Attenuation and stellar mass-related quantities 

By using optical-NIR flux densities, we can estimate the stel- 
lar mass of the sample galaxies. Here we discuss dust attenua- 
tion properties with respect to the stellar mass. There are vari- 
ous methods to estimate stellar m ass. In this work, w e used an 
SDSS-based method proposed by Yang et al.l d2007l) : these au- 
thors assume the Kroupa IMF ( Kroupal 1200 ll) . while we used 
the Salpeter IMF dSalpeterMl955l) to estimate the SFR of galax- 
ies. For consistency, we conve rt their stel l ar mas s estimates using 
the conversion factor given bv lBell et al.l d2003l) . As for the accu- 
racy of the stellar mass estimates, the relation has been calibrated 
from SDSS+2MASS data, and the error should be around 20 %. 
We also estimat ed the stellar mas s from Ks-band using the rela- 
tions derived bv lBell et al.l d2003l) . and the trends did not change 
significantly. Further exploration for the stellar mass estimation 
will be discussed elsewhere. 

Figure [16] shows the relation between stellar mass M* and 
dust attenuation in terms of Ljjr/Lfuv- At a first glance, the 
dependence of dust attenuation on stellar mass is very strong. 
Though we should be cautious on the selection effect on this 
plot again, the strong dependence might be partially physical. 
This implies that larger galaxies are more extinguished. Hence, 
dust attenuation is closely related to the physical size of galaxies. 

The presence of low mass and high Ljir/Lfuv galaxies 
~ 10-100 at M, ^ 10 7 — 1 9 M Q is worth mentioning here. 

esias-Paramo et al. d2006l) made essentially the same analy- 
sis by GALEX and IRAS. In their studies, such galaxies did not 
exist in their IR-selected samples (their Figure 12). If the flux 
measurement is secure, these galaxies deserve a close inspection 
to examine their nature. 

Next, we show the relation between SFR per unit stellar 
mass, SFR/M t , which is often referred to as the specific SFR 
(SSFR), and dust attenuation Ljir/Lfuv in Figure [17] It is in- 
deed striking that obviously there is no correlation between these 
quantities. 

First conclusion would be that there is no correlation be- 
tween global attenuation and SSFR. Generally, many authors 
showed tha t dust attenuation increases with SFR (see, e.g., 
Figure 7 of lBuat et all l2007al) . Then, if there is no link with 
the SSFR it implies that dust attenuation is linked to the total 
amount of SF scaled with galaxy size, because SFR was posi- 
tively correlated to M, seen above. It may be interpreted as fol- 
lows: the attenuation is not related to the current-to-past SFR ra- 
tio, or roughly, the SF history. Probably this could be explained 
by the short lifetime of dust grains in the ISM. 

We should mention that, however, esias-Paramo et al. 
d2006l) has shown that even for the FIR-selected galaxies there 



is some correlation between these two quantities (see their 
Figure 10b). This issue seems to require further examination. 



5. Conclusion 

In order to explore SEDs, star formation, and dust extinction 
properties of galaxies in the Local Universe, we have constructed 
a multiband galaxy sample based on the AKARI FIS All-Sky 
Survey and GALEX All-Sky Imaging Survey (AIS). We start 
from AKARI All-Sky Survey Bright Source Catalogue /3-1.1, 
and selected galaxies by matching the AKARI sources with 
those of the IRAS PSCz. Next, we have measured total GALEX 
FUV and NUV flux densities by a photometry software which 
we have developed specifically for this purpose. Then, we have 
matched this sample with SDSS and 2MASS galaxies to ob- 
tain the basic sample. The basic sample consists of 776 galax- 
ies. After removing objects with photometry contaminated by 
foregound sources (mainly in SDSS), we have defined the "se- 
cure sample" which contains 607 galaxies. Based on this galaxy 
sample, we have explored various properties of galaxies related 
to star formation and dust extinction. 

Summary and conclusions of this study are as follows: 

1. The sample galaxies have redshifts ;S 0.15, and their 90- 
fim luminosities range from 10 6 to 10 12 L , with a peak at 
1O 1O L . 

2. The SEDs display a very large variety, especially more than 
four orders of magnitude at M-FIR, but if we sort the sample 
by 90 pm, their average SED has a coherent trend: the more 
luminous at 90 /im, the redder the global SED becomes. 

3. The M, -NUV - r color-magnitude relation of our sample 
does not show a bimodality which is almost always expected 
in optically selected galaxy samples. The distribution is uni- 
modal, centered on the green valley between blue cloud and 
red sequence seen in optical surveys. 

4. We have established formulae to convert FIR luminosity 
from AKARI bands to the total infrared (IR) luminosity Ljir- 

5. The luminosity related to star formation activity (L$f) is 
dominated by the contribution of Ljir even if we take into 
account the FIR emission from dust heated by old stars. 

6. With these formulae, we calculated the star formation di- 
rectly visible with FUV and hidden by dust. At high star for- 
mation rate (SFR) (> 20 M yr~'), the fraction of directly 
visible SFR, SFRpuv, decreases. 

7. We estimated the ratio, Ljir/Lfuv, which is a direct measure 
of the FUV attenuation Apuv- The distribution of Ljir/Lfuv 
is con sistent with a prev ious result based on GALEX and 
IRAS dBuatet aUr2007al) . 

8. We also examined the L T ir/Lfuv-UV slope (FUV - NUV) 
relation. The majority of the sample has Ljir/Lfuv ratios 
5 to 10 times lower than e xpected from the local starburst 
relation dKong et al.L |2004|) . while some LIRGs and all the 
ULIRGs of this sample have higher Ljir/Lfuv ratios. This 
trend was already repor ted from a previous GALEX-IRAS 
study dBuat et al.Ll2005l) obtained by a much smaller sample, 
and we have confirmed their conclusion. 

9. By making use of stellar mass information derived from 
SDSS flux densities in this work), we have examined the dust 
attenuation properties in terms of stellar mass. We found that 
the attenuation indicator Ljir/Lfuv is correlated to stellar 
mass of galaxies, M„, but there is no correlation with spe- 
cific SFR (SSFR), SFR/M*. This may mean that Ljir/Lfuv 
is not linked to the SF history, but simply scales with the size 
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Fig. 15. Distribution of Ltir-^fuv luminosity ratio as a function of UV color FUV - NUV. Dots: galaxies with 60yum luminosity 
L(,o < 10 11 Lq, filled squares: galaxies with 10" < < 10 12 L (IR luminous galaxies: LIRGs), and filled circles: galaxies with 
10 11 < L(,q < 10 12 Lq (Ultraluminous IR l uminous galaxies: ULIRGs). Solid lines in each panel represent the relation for Local 
starbursts proposed by iMeurer et al.l dl999h . which was converted into the relation between the FUV-NUV color [eq. (fT9]ll. Left 
panel shows the raw distribution of Ltir-Lfuv luminosity ratio, while in Right panel functional fits are overplotted on the data. For 
galaxies with < 10 1 1 L Q , eq. ( TT9b is used, while for more luminous galaxies, simple linear relations are adopted because of poor 
fit of the function. Dotted line: a fit to galaxies with L( ,q < 10 11 L by eq. d20i i: Dot-dot-dot-dashed line: a fit to a sample of normal 
galaxies by eq. < f2Qb presented bv lBoissier et alj d2007l) : Dashed line: a linear fit to LIRGs; Dot-dashed line: a linear fit to ULIRGs. 



of galaxies. However, this is at odds with previous result of 
llglesias-Paramo et all d2006l) . 

This sample will serve as an important reference sample at 
z — for various further analysis or ongoing/future observational 
projects, like Herschel: for instance, this can be used to construct 
a set of SEDs for discussing higher-z observational strategy, or 
as a baseline test sample to investigate a m ethod of extractin g 
galaxies only from FIR flux information (e.g. JPollo et aUl20ld) . 

However, since our first sample is not complete in many 
senses, further analysis will be desired. We plan to construct a 
larger multiwavelength sample from the next release of AKARI 
FIS All-Sky Survey in the near future. 
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Fig. A.l. Luminosity vs. redshift diagrams for all the bands in the sample. 



